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Abstract 



Incompressible (magic) states that result from many-body effects in ver- 
tically coupled quantum dots submitted to strong magnetic fields such 
that only the lowest Landau level is relevant are studied within an exact 
diagonalization calculation for = 3, 5 and 6, electrons. We find that 
the sequences of total angular momentum M for which these incompress- 
ible states exist depend on the interplay between the inter-dot hopping 
parameter and the inter-dot distance d. For d of the order of the 
magnetic length and for all values of A^, we conclude that, in contrast 
to previous claims, these incompressible states appear at magic values 
of M which do not differ from those obtained for a single dot, namely 
M = N{N — l)/2 + j N where j is a positive integer number. For large 
inter-dot distance and simultaneously small inter-dot hopping parameter, 
new sequences of magic values of M are observed. These new sequences 
can be easily understood in terms of a transition regime towards a system 
of two decoupled single dots. However, important differences in the nature 
of the incompressible ground states are found with respect to those of a 
single dot. 
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I. Introduction 



Much effort has been devoted to understand the magic incompressible states (IS's) of 
two dimensional electronic nanostructures. This is due to the fact that they are closely 
related to the states that determine properties like superconductivity or the quantum 
Hall effect (QHE)^-^, which are striking examples of the non-trivial behavior that strongly 
interacting electronic systems may display^'^. Finite systems like quantum dots (QD's) 
provide simpler physical realizations of strongly interacting electronic systems where 
different models can be tested. When they are submitted to strong magnetic fields, the 
projection of the system to the lowest Landau level (LLL) becomes a good approximation 
which greatly simplifies theoretical studies in general, and, in particular, makes exact 
diagonalization calculations feasible. Much work has been done on single QD's in the 
LLL regime yielding a reasonable understanding of the nature of their IS's^'^. The search 
of IS with well defined properties which may produce fractional QHE experimentally 
observable, led to analyze double layered systems Double quantum dots (DQD's) 
in a vertical configuration submitted to strong magnetic fields provide a finite system in 
which the existence of IS is expected. However, the additional degree of freedom, together 
with the two new parameters, namely the distance between the dots and the tunneling 
strength, may give rise to new phenomenology. For instance, Yang et al.^^ suggest an 
experiment to test the quantum coherence of a special stable two-level system built in 
a DQD submitted to an adjustable interlayer bias voltage, which demonstrates suitable 
conditions for serving as quantum computing bits. Moreover, correlation effects can be 
experimentally detected in the far infrared range (FIR) using uniform electric fields with 
non- vanishing component along the vertical direction as the generalized Kohn theorem, 
under such condition does not apply^^. 

This paper is organized as follows. In Section II we describe the model used in 
our calculation and analyze the Hamiltonian of the system. In Section III, after the 
identification of the incompressible states of interacting electrons, we begin with a review 
of the results previously obtained for single dots and show next our main results for 
double dots, which cover a wide range of input parameters. Finally, in Section IV we 
compare our findings with previous results in the literature and draw our conclusions. 



II. The Hamiltonian 

We consider two identical two-dimensional quantum dots (in a vertical configura- 
tion) confined to the XY-plane by equal parabolic potentials and submitted to a strong 
magnetic field directed along an arbitrary direction. The Hamiltonian of the system 
reads, 

H = Ho + Ht + He-e (1) 

where Hq is the single-particle part which contains the kinetic contribution, the confining 
potential and the Zeeman term. We adjust the input parameters in such a way that 
Landau level mixing is negligible. Then, in second quantization formalism is given by 
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Ho^ aM + - AzS 



(2) 



where 

h 



«= 2(V^c+4c^o-^c) , (3) 



P^^yju;! + ^u;i , (4) 
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and 



IJ'sgB (5) 



ujq being the confining potential frequency, Uc the cyclotron frequency given by Uc = 
eB/m*c (m* is the effective electron mass, B the magnetic field and e and c the electron 
charge and the speed of light in vacuum respectively), Hb — eh/2mc the Bohr magneton 
and g the Lande factor (we will consider | g' |= 0.44 whenever the Zeeman term is 
included). M — J2iLi iTiia^.a^i is the total angular momentum and is the total number 
of electrons. a+. creates a single particle state and (Xj refers to the three indexes that 
characterize the single particle wave functions: angular momentum, spin and isospin 
( s or a associated to symmetric and antisymmetric combinations of wave functions 
concentrated in each dot: right and left). The tunneling term is given by 

Ht--Y^ (6) 

where At is the energy gap between the symmetric and antisymmetric states in the 
noninteracting system and X — Ns — Na is given by the balance between symmetric 
and antisymmetric states. Finally the two-body interaction part of the Hamiltonian is 
given by 

He-e = ^ E E Vm , V^] = {ij I V^^^ I kl) (7) 

^ ijkl A=l 

where the index A is used to distinguish between the three different possibilities: (i) 
V^^^ = when only one change of a single particle isospin takes place, (ii) V^"^^ = 
^{Vrr + Vri) wheu both isospins remain unchanged and (iii) V^^^ = ^(Kr — Vri) when 
both isospins are changed^^. Kr- and Vri are the intra and inter-dot Coulomb potentials 
respectively, which are given by 

Vrr = - (8) 

er 

and 
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d being the distance between the dots along the 2;-direction, r a 2-dimensional vector 
and € is the dielectric constant of the host semiconductor. We have assumed Dirac- 

delta distributions along the z-direction and have taken as a basis, Slater determinants 
built up from Foch-Darwin single particle wave functions projected on the LLL^. The 
diagonalization can be performed in separated subspaces characterized by three well 
defined quantum numbers: the total angular momentum M, the total spin S along the 

— * 

direction of the field B and the parity P related to the reflexion symmetry with respect 

to the plane midway between the dots ( P defined as P = (—1)"^/^ for even A^" and 
P = (— for odd N ). We wiU define the set (M, S, P) as a configuration. 

The eigenstates within each configuration are determined by H^-e + Ht alone, and 
the role of the constant term given by Hq is to shift the eigenenergies as a whole without 
changing their relative order. 



III. Incompressible states in the LLL. 



A. Single QD. 



Before studying the IS's in DQD's, we briefly review previous work on single QD's 
and its consequences. For a QD an IS with total energy E and characterized by (M, S) is 
identifled as the one which has the following singular property^^: the lowest excited state 
with quantum numbers (M+1, S) has energy E+a. That is to say, the energetically most 
favorable way to excite an IS increasing its total angular momentum by one unit is by 
moving the system as a whole, namely by increasing by one unit the angular momentum 
of the center of mass (CM) only and leaving the internal structure unchanged. This 
characteristic was nicely recognized analyzing the Coulomb contribution to the total 
energy of a full polarized QD as a function of M. A periodical arrangement of plateaux 
(steplike structure) in the otherwise decreasing curve signaled the values of the magic 
angular momenta^^. Furthermore, the variation of the magnetic fleld (or the conflning 
potential) did not drive the ground state (GS) through all neighboring values of M but 
through the sequence of magic values only^^. 

This scenario corresponds to the regime characterized by a flUing factor lower or equal 
to one, defined as^^ 

N(N-l) , , 

V = — 10 

2M ^ ' 

which involves the minimum possible value of the total angular momentum for a full 
polarized QD given by Mmin = N(N — l)/2 (the "compact state") and the angular 
momentum M of the magic state. Some care must be taken for low values of B for which 
the assumption of the LLL regime is not fulfilled. A suitable way to check this condition 
is by making sure that the energy of the highest single-particle occupied state is much 



smaller than a;+ = t( y + 4:^o+^c ) ; which is the energy gap between Landau levels 
for non-interacting electrons. The fiUing factor refers to the number of sublevels occupied 
within the LLL. There are two sublevels (spin up and down) in the case of a single QD 
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and four (two for spin and two for isospin) in the case of a DQD. In general, for regimes 
in which several sublevels are occupied, the fiUing factor of a QD is not well defined. GS's 
which are not related with IS are also possible under such multiple-sublevel occupancy. 

The sequence of magic filling factors depend on the number of electrons, for N ~ 3 
the values of u are z/ = 1, |, |, j, .. or for = 4 they are u — 1,^, |, |, .., in both cases 
related to the magic angular momentum given by 

M ^-N{N -1) + jN (11) 
2 

where j is a positive integer number. It turns out that the analysis of the Coulomb 
contribution to the total energy as a function of M gives exhaustive and precise infor- 
mation about the magic values of the angular momentum and hence about the magic 
filling factors. The magic values of M are the initial values of the plateaux. However, 
no information about the total spin of the IS's comes from the previous analysis. In the 
N — 3 case, for a QD, the sequence of GS's is always full polarized {Sz = 3/2) if the 
Zeeman term is included in the Hamiltonian (with | (? |= 0.44 ) or in contrast, oscilla- 
tions between = 3/2 and Sz = 1/2 were obtained if no Zeeman term is included in the 
calculation^^'^^. However, in the last case, the changes in spin and angular momentum 
do not appear simultaneously. 



B. DQD for d~ Is. 



For a DQD we have a richer parameter space to be explored as, in addition to 
the parameters of a single QD, and d also enter the Hamiltonian, which open new 
possibihties for IS's to exist. We will focus on the phase diagram (A^ / d) for standard 
values of the remaining input parameters. Due to the fact that Coulomb interaction and 
changes in parity are coupled processes in a DQD, we define the "interaction" energy as 
the Coulomb plus the tunneling contribution (C-l-T). 

For d Ib, where Ib is the magnetic length given by Ib = \j m*{ujl+Au^ I^^ P"^® 

Coulomb contribution to the total energy ( a = /3 = At = A^ = 0) as a function of M, 
is a decreasing function without plateaux as it is shown in the curves (a) in Fig.l and 2 
for = 5 and TV = 6 respectively (energies are given in units oi u = e^/els)- Figs.lA 
and 2 A correspond to parity P — 1 and Figs. IB and 2B to P = —1. All four cases 
refer to full polarized systems (S=5/2 for N=5 and S=3 for N=6) . For each value of 
M, the energy displayed is the lowest within the configuration (M, S, P) . The absence 
of plateaux can be understood as follows. Since A^ = 0, the number of electrons in each 
dot is a well defined number. Hence, in order to increase the total angular momentum 
by one unit, the angular momentum of either dot must be increased by one unit, which 
unavoidably increases the typical distance from the electrons of one dot to the ones of 
the other dot, and, therefore, decreases the inter-dot Coulomb energy. 

According to Figs. 1 and 2, it is necessary to include a sizeable tunneling contribution 
in order to obtain a sequence of plateaux, which, furthermore, only occur for P = — 1. 
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Indeed, from a series of calculations for N — 5 (not shown in Figs.l), which correspond to 
a variation of from to 2 meV by small steps, we see a number of plateaux gradually 
appearing as increases. We find that from At = 2 meV to 0.8 meV the sections from 
M = 10 to M = 11 and from M = 15 to M = 16 are exact plateaux. For A^ = 0.4 meV 
they are approximately fiat and for A^ = 0.2 meV they disappear. However in all cases 
the curves are abruptly decreasing before M — 10 and between magic values. That is to 
say, we do not find any extra value of magic M different from those given by Eq.(ll). 

In order to obtain information about the spin and parity of the IS's, we calculated the 
Coulomb plus tunneling contribution for all possible configurations. Fig. 3 (for N = 5) 
shows that the sequence of plateaux appear only when the system is fully polarized in 
spin and have parity P = —1 (similar results were obtained for = 6). Furthermore, 
although the parity P = —1 for N=5 can be obtained from different values of X, i.e., 
X — Ns — Na = 5, 1 or -3, the occupancy of the single-particle states for such 
incompressible GS turns out to be X = 5 only, namely the system is always full spin 
and isospin polarized. This suggests that the GS that are IS's will not present variations 
in 5* or P as i? increases. This last suggestion was confirmed, for Zeeman contribution 
different from zero, by an explicit calculation of the GS vs B, which turns out to be always 
full spin and isospin polarized. Fig.4 displays the total energy of the GS as a function of 
the magnetic field. The arrows point to the places where the angular momentum jumps 
from one magic value to the next one, leaving the spin and parity unchanged. In the 
inset we show Eos — (3N in order to compare with other publications which omit the 
A'^-dependent term. The nearly monotonous function of B is due to the fact that, in the 
absence of spin or isospin transitions, the interaction energy has a negligible infiuence in 
the plot and hence the evolution of the system is driven by the monotonous increasing 
term j3N which is much more important than the decreasing term aM, which would 
produce kinks at the transition points, as it is shown in the inset {Eos — I3N vs B). In 
brief, the full spin and isospin polarization appears to be a well defined attribute of these 
IS's that result from many-body effects. 

As shown in Fig. 3A the interaction energy appears to be degenerated at the magic 
values with respect to the three possible spin polarizations. Since the curves that belong 
to 5" = I and | have lower energy at the end of the plateaux, the final balance of energy 
depends critically on the relation between this difference of interaction energy, the kinetic 
and the Zeeman terms. That is to say, an IS that is the GS for a given value of B and 
ujq will remain as GS as B increases or ujq decreases only if 

Ec+t{M, ^,1) - Ec+t{M + 1,^,1) < a + Az (12) 

or otherwise the new GS will be a compressible not fully spin polarized state at M + 1. 
Hence, as B or ujq change, the GS can be driven into compressible zones in contrast with 
the results obtained for single QD's. 

The single particle occupancies of the m-values for the first three IS's for N = 5, 
calculated at S = 4, 7 and 9 T respectively are shown in Fig. 5. The first GS for M = 10 
(i/ = 1) is the compact full polarized state which belongs to a one dimensional subspace 
and, as a consequence, no correlation is involved as one Slater determinant produces the 
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exact solution. Moreover, the density is a "dome" shape circular symmetric distribution 
without any structure. As B grows, the angular momentum changes from M to M + A^, 

all the electrons jump together moving away from the origin and forming a ring. The 
dimension of the GS subspace increases and as the relative weights of the different Slater 
determinants within the expansion of the GS become significant for different elements 
of the bases, namely, the correlation becomes important. 

C. DQD: Phase diagram At/d. 

So far we have explored the situation d ~ Ib- In order to get the complete scenario 
of IS's in a DQD, we have also investigated in detail the remaining regions of the phase 
diagram (At/d) 

In Figs. 6-9 we follow, for N — 3, the variation of the interaction energy vs M {E{C + 
T) /M) as d and At change. 

Prom A to B (Fig. 6), as it was just discussed for the case d ^ Ib, the plateaux 
emerge as grows from zero until they are well defined at A^ = 0.11 -u, that is to say, 
for values of M given by Eq.(ll) the system evolves from compressible to incompressible 
states. For d close to zero (at the left of point B) and Aj large, all the electrons are in 
the symmetric state. As a consequence, the interaction energy of a single QD can be 
reproduced with high accuracy by the addition of the constant contribution A(X/2 to 
the energy of the DQD in this region. 

From A to D (Fig. 7), tunneling between the two dots is not allowed. Starting from 
a curve without plateaux for small distances {d ~ lOA), we move across the transition 
regime with a gradual formation of new plateaux at M = Mr + where Mr and Mj, 
are the magic numbers of single QD's. We come close to the point D at d = lOOOA 
which show the features of two decoupled dots with = 1 (with no contribution to 
the Coulomb term) and N = 2 (with magic numbers M = 1, 3, 5, 7..) respectively. This 
regime in which tunneling is forbidden has been previously studied for double layers 
and special attention has been devoted to the v = \ case^'^'^^. For a double layer 
the incompressible state v — ISs observed for values of d about the magnetic length^. 
Furthermore, as d increases, the state exhibits a phase transition to a compressible one. 
The difference between the two cases comes from the fact that, as it was previously 
discussed, when Af ~ Mr and are well defined quantum numbers, the increase by 
one unit of M means the increase of Mb or Mb (but not both) changing, in a DQD, the 
relative position of the electrons in each dot and so decreasing the Coulomb interaction 
which prevents the formation of a plateau. However, this is not the case for a double layer 
in which the shift of charge due to the change of angular momentum does not change 
the relative inter-layer distribution of charge, allowing for the appearance of plateaus. 
For large d the two layers decouple and hence one would expect (for total v = \) two 
V — 1/2 IS's. However, since fractional QHE oi v — 1/2 is not observable for a single 
layer, these states were not identified in Ref.(^). 

From B to C (Fig.8), although the distance d grows, the sequence of magic numbers 
typical of a DQD does not disappear due to the relative large value of A^ ( A^ = 0.229 u) . 
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An exceptional case appears for d > 500A at M = 1 (for N — 5 the analog exceptional 
case appears at M = 4). It turns out to be the only IS within the LLL regime which 
does not fulfill the general rule of being full spin and isospin polarized. The subspace 
associated to the appropriate configuration, i.e., (M = 1, S* = 3/2, P = 1) is one dimen- 
sional and the only Slater determinant in the bases has one electron in the symmetric 
state and two electrons in the antisymmetric state (or X — —1). It is the only IS with 
no single dot analog. Our interpretation of the fact that the M — 1 magic value appears 
only for relative large distances is as follows: for large tunneling and small d, the system 
is closer to a single dot with = 3 than to a DQD of the same number of electrons (as 
discussed before). Thus M = 1 can only appear when the Coulomb inter-dot interaction 
weakens related to and DQD properties different from those of a single QD may 
arise. Notice however, that for the M — 1 state to be a GS such a low magnetic field (or 
large confining potential) is required that the LLL regime assumption would not apply 
anymore. Finally, even for values of d as large as 1500 A (being Ib = 65^4), we did not 
find the transition from DQD to two decoupled single dots. 

From D to C (Fig. 9), the tunneling increases and the system of two decoupled QD's 
with a period of two typical of the N = 2 single dot evolves into a DQD reproducing 
the period of three typical of a = 3 DQD. During the transition, there is a narrow 
interval of values of At for which the E{C + T) ys M curve has no plateaux (except 
for the M = 1 case). That is to say, an initially incompressible GS would evolve into a 
compressible state and again into a IS as At increases. This evolution takes place as the 
system changes from two decoupled single dots to a DQD. 

Compressible regions have been obtained before by Rontani et al.^^ for DQD with 
finite width. They consider the evolution of the GS of the system of A^ = 6 as d 
increases for At exponentially decreasing with d, which is equivalent to the evolution 
along a trajectory from B to D in our phase diagram. They obtain a small zone of 
compressibility in the middle, related to the transition from a regime where the system 
behaves as a unique coherent system to a regime of well separated QD's. We observe 
the same qualitative behavior along the B-C trajectory (which is different to theirs) 
although we do not obtain the same magic values. 

The transition from DQD to two decoupled single QD's as At decreases has been 
observed before by Peeters et al^^ by means of a current spin density functional calcula- 
tion. 

We have also studied the D to C evolution for N=5 which shows the same quahtative 
behavior. As D is approached the structure of plateaux can be understood in terms of 
the IS of two decoupled N = 2 and A^ = 3 single QD's, although the analysis is much 
more intricate than for the N=3 case. 



IV. Discussion 

We have investigated in detail the existence of IS's that result from the Coulomb 
many-body effects in a DQD for the entire phase diagram (At / d). 
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An important point in our analysis is the criterium used to identify IS's. First we 
want to emphasize that in contrast to the case of a single layer for which the integer 

QHE is associated to gaps of single-particic origin and the fractional QHE is associated to 
gaps involving many-body effects, for double layered systems (and thus in accordance for 
DQD), single-particle as well as many-body regimes can be related to QHE at the same 
filling factor by the tuning of appropriate sample parameters^. The IS's we are interested 
in are those associated with e-e interaction (coupled with tunneling) and thus signaled 
somehow in the variation of the interaction energy with M. As discussed in the previous 
section, we define the interaction energy as the Coulomb plus tunneling contributions, 
and require IS's to preserve the interaction energy when the angular momentum M is 
increased by one unit. We want to stress that this is not equivalent to identifying magic 
M from the kinks of the lowest energies of each configuration as a function of M or from 
the kinks of the variation of the absolute GS energy as a function of B, as it has been 
used in the literature^^''^^ to identify correlated IS's. 

Our criterium is equivalent to the one used by Laughlin in Rcf.(^) for the single 
layer, as we discuss below. The Hamiltonian is separable into the CM and the rela- 
tive coordinates and as a consequence, the total angular momentum can be analyzed 
as M = McM + Mr el and the energy as Etot — Ecu + Einternai- For three two dimen- 
sional electrons, Laughlin obtains that the internal energy as a function of the relative 
angular momentum has downward cusps at special (magic) values (M = 3,6,9,12,..). 
These magic values appear to be related to incompressibility: the area of the system 
defined as the area of the triangle determined by the correlated positions of the elec- 
trons within these states changes discontinuously as pressure is applied. At the down- 
ward cusps, EinternaiiMrei) < Einternai{Mrei + 1) for the lowcst energy statcs of cach 
configuration. They arc the only states for which the increase of Mrd by one unit 
requires a positive amount of internal energy. In order to show that our criterium^^ 
is equivalent to Laughlin's^, notice first that Ec+t{M) only depends on M^ei and 
Mj-ei < M . Since Ec+t{M) is defined as the minimal energy among those of the 
states with total angular momentum M, it implies it is the minimal energy among all 
states with relative angular momentum Mrei < M. Hence, given M and Ec+t{M — 1), 
Ec+t{M-1) ^ Ec+t{M) implies Ec+t{M-1) > Ec+t{M) and furthermore M = M^ei. 
Since Ec+t reduces to Einternai for a single layer, if Ec+t{M — 1) > Ec+t{M) then 
EinternaiiMrei " 1) > EinternaiiMrei) for Mrei = M. Namely, negative slopcs in our plots 
imply negative slopes in Laughlin's. If, on the contrary, Ec+riM — 1) = Ec+riM), 
then Mrei has not changed and M contains at least one unit of CM angular momen- 
tum. Being Ec+riM) the minimal energy with total angular momentum M, it implies 
that any state with Mrei = M has larger energy than Ec+riM — 1). If, in addition, 
Ec+riM — 2) > Ec+riM — 1), as it is always the case in our plots, then M — 1 con- 
tains only relative angular momentum, as shown above. Then, for the single layer, 

Ec+riM - 1) = Ec+riM) implies EinternaiiMrei - 1) < EinternaiiMrei) for Mrei = M. 

Namely plateaux in our plots imply positive slopes in Laughlin's, which concludes our 
proof. On the other hand, there is no guarantee that downward cusps in the curve 
Etot vs M are related to plateaux in the curve Ec+t vs M. It is enough to have a > 
Ec+riM) — Ec+riM -|- 1) at a non-magic M to obtain there a downward peak in Etot 
vs M which is not related to a plateau in Ec+r vs M. 
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For d Ib the fact that the magic M follow Eq.(ll) is in conflict with the claim 
made in Ref.^^'^^ that extra magic values for M (depending on the value on the tunneling 
strength) exist in this regime. The authors of Rcf.(^^'^^) identified IS's with downward 
cusps of the total energy as a function of M, i.e. the interaction energy (our curve) with 
the addition of the single-particle contribution. In Fig. 10 we show the two possibilities 
for = 5. It is clear from the upper curve that some downward cusps, which would be 
identified as IS's by the criterium of Ref.^^'^^, do not actually correspond to IS's in our 
criterium. 

In order to make sure that the discrepancies with Ref.^^'^^ arc only due to the different 
criteria to identify IS's, we have reproduced their results (sec Fig. 11). To be more precise, 
we performed the calculation for = 3 and the same input parameters as those used in 
Ref.^^ { N = 3, B = d = 200 A, fiwo = SmeV and At = 0.2 meV ). In Fig.llA 

the interaction energy contribution versus M is shown. Due to the low value of the 
tunneling contribution, the plateaux that will appear, for larger values of A^, at M= 3, 
6 and 9 are still not visible and the only ones that already appear are M= 12 and 15. If 
the kinetic contribution is added, it comes out that the GS is at M= 5 (see Fig.llB) at 
the lowest downward cusp of the total energy in accordance with Ref.^^. 

A word of caution should be given here as a number of relevant papers exists in 
the literature^^"^^ which use the term magic angular momentum to denote the angular 
momentum M which displays downward cusps in the curve Etot{M). The corresponding 
states enjoy enhance stability and have been the subject of intensive studies. In particular 
Refs. (^^'^*') provide characterizations of these states ranging from small values of M, 
where the fractional quantum Hall regime is sometimes identified, to large values of M, 
where strip- like structures and Wigner molecules seem to appear (see Ref. (^^) for a 
review and Ref. (^^) for related work on layers). However, only a subset of these states 
fulfills our criterium of incompressibility and the magic M displayed in (11) correspond 
to this subset only. 

Our conclusions can be summarized as follows: 

(i) The downward cusps obtained by Laughlin in Ref.(^) turn out to be equivalent to 
the plateaux of the curve Ec+t vs M. 

(ii) All the incompressible states are full spin and isospin polarized (except the M = 1 
case for A^ = 3). Since a single QD full spin polarized and a DQD full spin ans isospin 
polarized are systems with no extra degrees of freedom aside from angular momentum, 
we expect a similar behavior for the electronic distance quantization as that obtained in 
Ref(i). 

(iii) An exceptional incompressible state was found iov d/ls > 8 at M = 1 and 
X = -1 for A^ = 3 and M = 4 and X = 1 for A^ = 5. This is the only one that is not 
full isospin polarized and docs not have a single QD analog. However, for it to be a GS 
values of the input parameters that do not fulfill the assumption of the LLL regime are 
required. 

(iv) For d ~ Ib, it is not possible to obtain IS's if the tunneling is small. 
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(v) For small A^, as d grows, the DQD evolves into two decoupled single QD's. New 
magic values appear which correspond to the addition of magic numbers of two decoupled 
QD's. 

(vi) An IS at M will remain as the GS under changes of B or/and loq if the condition 

Ec+t{M) - Ec+t{M + 1) < a + Az (13) 

is fulfilled. In general, however, the variation of S, cuq or can drive the GS from 
incompressible into compressible zones of the phase diagram. This behavior differs from 
the well known properties of a single QD for which the variation of the GS as B or ujq 
change is driven through incompressible states only, skiping all non-magic values of M. 

(vii) Whenever the GS is full spin and isospin polarized, it is an IS. In other words, 
the IS's are the only possible GS for u lower that one. However for i/ > 1 other GS are 
possible. 

Before closing, let us briefly elaborate on the last point. Notice that the following 
situations are also possible: (a) GS's with M that fulfills Eq. (11) and arc not IS's due 
to the fact that from M to M + 1 there is not a plateau. This condition can be obtained 
for very low values of A^, for example, for iV = 5, M = 10, B = 5T, Uujq = 2.6 meV, 
At = 0.2 meV and d = 20 A. (b) Ground 

states with M not given by Eq. (11) for which the system is not full spin or isospin 
polarized, (or z/ > 1 ) and E[C + T) has not a constant evolution from M to M + 1. 
This is the case for example for = 5, M = 13 ( -B = 6T, Hcuq = 2.6 meV d = 20A and 
At = 5.86 • 10~^meV). In the last case, S = N/2 and P = -1, however X ^ N. 

Notice also that the previously discussed states are not the only possible GS's within 
the LLL regime. For instance, if the confining potential is strong enough, other types of 
GS's are possible like the ferromagnetic, canted and symmetric states (all of them with 
u — 2) first studied in double layers^^ and latter recognized in DQD^^. 

Finally, let us note that the correlation plays an increasingly important role as the 
magnetic field grows up. The interaction energy and correlation effects can be experimen- 
tally tested by uniform electric fields with non-vanishing component along the z-direction 
due to the fact that under this condition, the Kohn theorem does not apply^^ and the 
FIR spectroscopy becomes sensitive to the internal structure. 

We gratefully acknowledge C. Tejedor and L. Martin-Moreno for the code used for 
the Hamiltonian diagonalization. This work has been performed under Grants No. 
BFM2002-01868 from DGESIC (Spain), No. FPA2001-3598 from MCyT and Feder 
(Spain) , and No. 2001GR-0064 and No. 2001SGR-00065 from Generahtat de Catalunya. 
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FIGURES 
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FIG. 1. A: Coulomb plus tunneling contribution to the total energy as a funtion of M for 
N = 5, S = N/2 and parity P = I, for several values of the tunneling gap: (a) = 0, (b) 
At = 2.2meV and (c) At = llmeV. B: The same as A for P = —1. The triangles point to 
the begining of the plateaux. We have taken B = 5T, d = 20A, and Hloq = 2.6 meV. 
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FIG. 3. A: Coulomb plus tunneling contribution to the total energy as a function of M for 
N = 5 and parity P = —1 for all the possible values of the spin S. B: The same as A for P = 1. 
We have taken B = 5T, d = 20l, Tnoo = 2.6meV and Aj = 2meV. 
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FIG. 4. A: Ground state energy versus B for N = 5. The triangles point to the transitions 
between magic angular momenta. We use the same values for d, ujq and At as in Fig.3. Inset: 
Egs — as a function of B. 
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FIG. 5. Total occupancy of the single-particle angular momentum states m for A?^ = 5. The 
values of B considered correspond to the magic values M = 10, 15 and 20 respectively. We 
have taken d = 20A , Huq = 2.6 meV and At = 6meV. 
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FIG. 6. Coulomb plus tunneling contribution versus M for iV = 3. Prom A to B the 
tunneling gap is: 0.0, 0.018, 0.106 and 0.229u respectively. Inset: phase diagram used. We 
have taken huo = 2.6 meF, B = 15T, d = lOl, S = 3/2 and P = 1 {Ib = 65l). 
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FIG. 10. Coulomb plus tunneling contribution (lower curve) and total energy (upper curve) 
as a function of M ioi N = 5. The black triangles point to the actual magic values M whereas 
the white triangles point to cusps which could be mistaken by them. The values of B, uo, d 
and A( are the same as in Fig.3 
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FIG. 11. A: Coulomb plus tunneling contribution versus M for A/" = 3. B: The same as in 
A for the total energy. The input parameters are given in the text. The black triangles point to 
the actual magic values M whereas the white triangles point to cusps which could be mistaken 
by them. 
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